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character ist ic iiietliQcJ ^artificial viscosity and yithom's rule * 


Trie second chapter is devoted to unify ' study of 
fuagnstcradiat ivs shock i^iave propaorat ion in conducting plasiia ^ 
taking line eKplosicn in a gas cloud and a detail study Is siads 
for the propagation of spherical and cylindrical shock, wave », 
tMunisrlcaii integration of differential equations of laotiari is dons 
on £u E « C i 0‘90 , cQ.fiipu ter at Iv I mT Kanpur by ws 1 1 Khoowri R « K » S « S 
Pr.ograrriine*. 


Cliaptar third. ' deals with -seif siniiar inagnetcgasclynaiiiic 
cylindrical shock waves with transverse and aoisutiial aKiai 
ffiagiietic field « The inotion of gas is assumed to 'tie, adiabatic and 













total energy of the wave remain cosntant . In this chapter 
instataneoas release of energy along a line in a gas cloud is 
asBused . A coffipar^at i ve study .is made between cylindrical shock 
wavs with transverse and aKiai magnetic field. J'iu.ffler icai 
intscration of differential equation's is carried out on DEU lu90 
co, 7 .puter system at I.I.T Kanpur by R.K.6.S progromme . 


in fourth chapter , an attempt is sfiads for the study of 
spharicaily syffimetricai str-ong discontiriuitiss with incrsaaing 
snsr-QV in generaiioed Roche model and numerical solution are also 
obtained when the radiation heat fiUK is more i moor- 1 ant then the 
radiation pressure and radiation energy . The effect of magnetic 
field has also been taken in to accout . 


Fifth chapter consists of analysis of seif similar' motion in 
the theory of stellar ss-iplosion, taking Newtonian gravitation into 
account, A thorough analysis of self similar equations of gas 
dynamics, under the effect of magnetic field is developed. 
A number of new solutions has been obtained in which radial 
oscillation of gas occur after the shock wave oasses -It is 
supposed tnat: star is in equilibrium state , 


efcudv of acf exact solutions of 


Sixth chao-tsr deals with ths 
r.oi-mai shock wave of varibais strength advancing into a region of 
variable density taking pressur's as constant and density vanes 

c-ccT.:' d I no to B.B poi4sr *. 


In seventh ciiapter.an attempt is made for the anaiyticai 
solution of cylindrical shock wave in a rotating gas with 
azimuthal magnetic field , taking density as constant ahead the 
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< i > SHOCK yAVES : - 

When disturbances of finite ampi i tude are 
oropaQated in perfect fluid or gases tfiat is those with 

no viscosity or heat conduct ivi tv ,, discont intii t ies in pressLire 
and veiacitv to the meditini may occur « These are called shock 
waves or shock . The reason for their deveiopmerit may be seen in 

the case- of one dimensional flow such as occurs in a tube of 

unitGrn cross section when a disturbance caused by the motion of 
a piston at one end, is propagated down the tubs* The principles of 
conservation of momentufn and energy is still apply, across , 

these plans discantinui t i.es«' These disturbances propagated 

through out the fluid as a wavs motion and with the speed of 

sound relative to the fluid without sufferiog any destoration* 
it may be pointed out that shock wave is a surface of 

discont iritii ty , ptilss like in nature and is some times more 
appropriattsiy called shock front C C 1 3, £23, £33» C43) . In general 
discont iriLiity surfaces are of two type •*- contact stirfacs and 
shock front A contact surface is a surface separating 

two parts of the medium with out any flow of a gas 

tfiroiigh the surface * A shock front is a discont iraiiity which 
is crossed by gas. 

The velocity Cu3 of the piston at any time is cofumunicated to 
the gas and propagates down the gas with velocity of soiirid Cc3 
relative to the gas- Hence relative to the tube the velocity is 
propagated with the, speed i 'c ^ u ) in the direction away 


i 



from the piston. if the piston moves to compress ths I 

i 


pas in the tube that is , to decrease the volume of the 
pas, then the velocity of the pas irsereases since the 
temperature and density increases. Thus the greater the velocity, 
the greater the speed with which it is propagated if this process 
continued indefinitely, lesser velocities would be over taken by 
greater ones and we would have two values of the velocity st a 
given place in the pas which is impossible. 

However, before this occurs there will be a time at which the 
velocity profile of the pas, that is the velocity distance curve, 
has a vertical shape. In this case, the differential equation 
poverninp the motion break down and the basis for the statement 
that the fluid velocity 1 c + u ), is no longer true - 

Before formulatinp these laws for a perfect fluid it may be 
pertinent to point out that it is consequence of the non - linear 
character of the equations . poverninp the oropaption of 
finite disturbance ; that is in part responsible for the 
prediction of the formation of the discontinui ties. Thus in 
the example mentioned above it is essential to the arpument that 
the local veiocit.yf of sound ic) be preater at points where the 
velocity of the gas is larger than velocity of sound. This is the 
case for a compressive motion of the piston because the 
disturbance which have passed over the pas have changed its 
character (when heated it) and hence subsequent disturbance are 
travelling in a defferent medium than the original one. 
Taking in to account . the change in the medium produced 
by one part of the phenomenon in the discussion above , 
subsequent part is accomplished mathemat icai ly by the non 


1 inear- 


terms in the equation of motion. 



me non 


. inear cnaraccsr aione is na-c enouian to cause 


discontinuities for it, in the example given above, the piston motion 
were such as to increase the volume of the gas in the tube, that 


is, if a rarefraction ^ave propagated down the tube then the local 


ve loci tv 


sound would decrease with increasing fluid velocity 


and even if a discontinuity were originally present it would 
disappear with time. The fact that only compression shocks are 
found in media which behave appraximateiy like ideal gases is in 
aqreeiment with the second law of thermodynamics : iWaiker and Taub 


£53) 


true of course, that a shock wave is not a discontinuity 


comp a t 


field. 


ickness acr 


perties change continuously. If this thic 
th some appropriate macroscopic dimensio 
as the radius of the curvature 
hysical relationship may be obtained b 


the discontinui tv as strict- The assump 


V thickness is small is a fundamental 


as applied to a shock wave, refers to 


opert: 


of the fluid within the sma 


if the discontinuity C63.If th 


ich the 


is small 


he f 1 ow 


curved 


nal vsii 


hat the 


he term 


aiue of 


finite 


ynamic 


a cnaraci:er is'cic time 


must elapse before squiiiburlum can be approximately 
reestablished. If the physical and chemical changes occuring in 
the discontinuity are sufficiently slow, so that the thickness of 


the discontinuity 


large compared with the characteristic 


d istartce, the concept of the thermodynamic quasi — equilibrium may 










Due to the high temperature there is of interest to consider 
the effects of thermal radiation in gasdynasjic. Sach Ciil was 
first person who derived shock conditions taking into account 
radiation pressure and radiation energy . Later, 6uess and Sen Ei23 
considered the effect of radiative transfer. Marshak Ei33 took 
into account radiation flUK also and he rederived the Hugnoit 
shock conditions and sifniiarity solutions . Elliot 11143 and Wang 
Ci53 considered the similarity solutions of Taylor's expansion 
probiem and Piston problem taking an account of Radiation flux. 
Bhatnagar and Sachdev tISl studied of the oropagation of 

isotherasal shock, thev used the Witham's rule to study the 
variations in the shock strength and shock velocity, after using 
pressure energy and fiuK, of radiation. The purpose of oresent 
work is to make thorough study of such radiative 
magneto— hvdrodynamic shocks by formulating suitable def ferential 
equations and reducing them to the form as to make them amenable 
for numerical integration. 

2. SPHERICAL AND CYLINDRICAL SHOCK 

Consider the propagation of a shock wave, through a perfect 
gas, of great intensity resulting from a strong explosion , ie from 
the instantaneous release of a large quantity of energy. 

When the energy is suddenly released, in an infinitely 

concentrated form and distribution of density oressore etc, 
depends only on the distance from some point then this is a case 
of spherical shock. 

When the energy is suddenly released along a line and 

distribution of all quantities is homogeneous in some direction 
and has complete axial symmetry about that direct ion, then this is 


e 





a case of cylindrical shock. 


EQUATiOMS OF MOTION AND JUMP CONDITION IN MAGNET! 


HYDRODYNAMICS 


eisctricai iv 


conductino fluids in the presence of 


maonetic field discont inui ty, ie. Shock wave in flow variable can 
exist "5.. The study of mapnetohydrodynasnic shock waves was bepun 
in 1950 with the paper of F.de Hoffmann and Tayiar LIOJ . 


When electric current induced in the fluid, then their flow in 
the ffiaonetic field produces mechanical forces which modifv the 
tnotion ti?!, maonetopas dynamic owes its peculiar interest and 
difficulty to this interaction between the field and the 

motion. Thus the equations of mapnetopasdynamic are the 

ordinary and slectromaanetic equations . Ne only take into 
account the interaction between the motion and the mapnetic 
field , we have ipnored the Maxwell's displacement currents. 
Throuph out the thesis we takeiu)— i because mapnetic permeability 
lu ) differs only siiphtly from unity which is unimoartant, We 
also assume that the dissioative mechanism such as 

viscosity, thermal conductivity and electrical resistant are absent, 


in th 


Since the problems dealt in this thesis relate to 
‘to gas dynamic or mapneto— radiative shocks . We refer 

:is article the relevant flow and fluid equation C1S3 
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where 


D 

Dt 


^ + y 

at 3>r 


and u,P,p,h,r and t are the velocity , pressure , density 
magneteic filed transverse to the flow , radial distance and 
tirae respectively , t is the ratio of specific heat at constant 
voluffie and constant pressure . -J=0. 1 and 2 corresponding to 
plane . cylindrical and spherieal respectively and V -U for plane 
and V =1 corresponds to cylindrical and spherical both . 


in' the presence of ®agnetic field, the relation connecting the 
Flow and field quantities on the two sides of the shock surface are 
m follows C 19, 203, where the velocity in front the shock wave is 


zero 


h (u-u > = h u 
2 2 1 


<u— u ) 


p u ? 
1 


1 




8 


p + 




Cii-u I = P 
2 1 


li +' p u 
1 1 


TP ii 

1 2 2 . 2 

— iLI““0 ) + 

2 2 CT“”i)p p 


TP tl 

. 2 ' !■ i 

li + + — “ 

CT-ii p f 


where sufiK i and 2 correspond to the value of the 
quantities iust ahead and just behind the shock surface and u is 
the shock velocity. 


(4) RAD I AT SON PHENOMENON 


At very high tsfiiperat-ure, the gases become ionized and 
radiation can be considered as a continuous emission of energy in 
the form of electromagnetic wave which propagated in the medium 
with the speed of light. This energy is called radiant or thermal 
energy radiation . Where as according to quantum theory the 
radiant energy emitted or absorbed is not continuous , 
periTiltting ail possible values . as demanded by the wave 
theroy , but in a discrete quantified form , as integral 
■ multi piss of an eiamentry quantum of energy , photon or 
light quanta .The amount of energy in each quanta being given by 
the product hO , where h is plank's constant and 9 is the 
frequency of the radiation. 


Thus the quantum theory proposes the particles caracter istic 
of radiation , while ciassicial theory the wave charactrist ics , 
both being required to understand the complex behaviour of 
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• . . .. 


radiation . A complete study of such a hiph temperature flow of a 
pas should consist of the study of pas dynamic field , the 
electromaonet ic field and thermal radiation flu;-; simultaneously . 
The anlvsis of simultaneously effects and magnetic field forms 
the Bub-iect matter of radiation magnetohydrodynamics - This 
theory of thermal radiation can be applied to understand the 
processes which take place in stellar media , to eKplain the 
observed luminosity of stars and nuclear sKpiosions and also to 
high temperature flow. 

Radiative transfer and radiative heat eKchange have an 
influence on both the star and the motion of the fluid .This 
influence is caused by the fact that fluid loses or pains energy 
by emitting or absorbing heat . On the flow field of the gas 
there are three radiation effects e;-;pressed in terms of radiation 
pressure , radiation energy and radiation fluK C21 ,22 ,233 - 

<a) RAD i AMT PRESSURE 


Bv the theory of electrodynamic the pressure of a radiation 
field is equal to one third of radiant energy is 


1 14 

P = — E = — a T ; 

R 3 r -3 R 


where T is temperature and aB is Stefan - Boltz mann constant - 
This is the only component of the radiation which differs fro® 
zero . 


( b > RAD i AT I ON ENERGY ( E, ) 


The radiant energy density E^. per unit mass of the fluid is 


given by 
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4 


a T 
R 

E 

R ^ 


where p is density of fluid 

ic) RADiftTI0t4 FLUX 

The net afRount of radiant sneroy passino throuph the sur f st-s 
per unit area per unit tiais is called the radiant fluK t-tU'Ouph 
the surface is aives by 


C 


R 4 


a T 
R 


or f = D prad E 
R R R 


CL 


where D" 


is the Rosssiand diffusion coefficient of 
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2 

h 

h 2 p 

where h is "fehs 0139061; ic ■field 

The ftsodaroen-tsi equation of radiation (aaor.stOQasdvr.amics for 
one d isnens iona i t 1 cik arS' C 24 3 


+ i = o 

Dt ^ hr r 


Du 

+ 

1 h Bh 

= 0 5 

Dt 

^ Sr 0 3 r 

pr 

Dh 

a u '^hu. 

« 

Dt 

^ r r 


DE 

Dt 

D 1 ^ 

Dt p j 

S J 

< Fr ; 

Ur 


D d a 

Dt 2>t. dr 

) ; 

where 

P = p +p +p 5 

SIS r h 
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T 



and F is the radiation heat FIuk. 

where suffixes fflt B.t h attached to the quantities to expressions 
far fiiaterial radiation and maonstic teras resoectivsly . 

5. SrHILARITY PRINCIPLE 

The fluid is said to be one-dimensional when all its 
properties depend on only one oeosnatric coordinate and on the 
time. The sohericai, cylindrical and plane waves produce one 
dimensional motion. The methods of dimensional analysis can be 
used to find exact solutions of certain problems of one 
dimensional unsteady motion of a compressible fluid. In the 
Eculerian appr^.och the basic physical variables are the velocity, u 
the density ^ and pressure P. The charactsrist ic parameters are 
the linear coordinate Cr3,the time t and the constants that enter 
into the equctt ions, the boundary and the initial conditions of the 
problem. Since the dimensions of the Quantities ^ and P contain 
the mass, at least one constant a, the dimension of which also 
contain the mass must be a characterist ic parameter. Hence, as in 
£241 we can assume with out any loss of generality 




o one* 
i-iB can 


irite thi 


.density and pressure 


i^iiiers 


V,6 arsd P are arbitral 


theref are, can depend only on one-dimensionai c 


and other d 


irasseter of the problems. In 


functions of tma non-dimensional variables. But 
characteristic parameter 'b' can be introduce 


irid&psndent of those o- 


The riLiifiber of iiidi: 


wiiicti can be formed by coffibinino a and fe i 

Since the dimension of the constant "a^ contair 
choose the constant 'b' in a manner. sucn thai 


not contain thg pfa 
dim b = L T 


Tlie s i nci i e iiari— d i insriS i ona 1 i ndependsn t var 


14111 oe 


can be replaced by m^O 


'Of! 










if and 


is o r'opQ^’^'^ ^ 


d p deDsnd only on the tifie and the veiocitv u 


siutions dependina 


on the independsni 


fariab les 


sv obtain a number of arbitrary constants 


The above aroument shows that when the characteristic 

parameters inciude two constant with the independent dimentions 


in addition to 


and t 


bhe partial def fsrieant ia i equation 


jatisfied bvthe velocity, density and pressure in one dimension 
unstedv motion as a compressib ie fluid can be replaced bv 


ordinary differential epuatians for V,6 and P. The solutions ot 
these ordinary differential equations can some times be obtained 
eKsctly in closed form £251. For violent spherical sKpiosions in a 
uniform atmosphere at rest , and in other cases, approKimated 
by using numerical integration £261 such motion are called seif 
similar £271. The idealised problems of a strong eKolosion ^.n 
homogeneous atmosphere represents a, typical eKamp le . of 

self-similar flow, in which the flow variables changes with times in 
such a manner that their distributions with respect to coordinate 
variable always remain similar in time. The self similar problem 
of a strong sKplosion was formulated and solved by Sedov £281 using 
a brilliant method, which employed the energy integral. Sedov, 
succeded in finding an exact analytic solution to the equation of 
self-simiiiar motion . The same problem was also considered by 
Tayic.r£2:63 who formulated the equation for the problem obtained 
nj'nericsi but not analytic solutions. 


eauat ions 
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linear in first derivatives .of the dependent variables, but the 

coefficient fsay be functions of the dependent variable when these 

epustiorss describe wave oiotion. s Qood understand ino of Hiany of 

the issues can be developed fross the study of plane wave=t« 

AccordinGiy,we start with the case of two independent variables 

The two variables often the time and space variable so we 

denote them by t and x , and use corrsspondino terminoicoy , but the 

discussion aoplies to any variable system if the dependent 

variafeiss are L* i = l,2 n, the penerai Quasi-linsar first 

i 

order system is 


3 u a u 

p. + 

id 9 1 i i 9 X 


+ b 


0 


(i. 1) 


where the matrices ft, a and vector b may be functions of v v 

,i n 

as well as ;< and t. 

In 'perteral ,any one of the equations has different combination 

of au /'at and /^>t far u . That is it couples information about 
i j J 

the rata of chanpe of the different u. in different directions . 


One can not deduce information about the increments of ail the u 

•j 

for a step in any sinqis direction but we are at liberty to 
manipulate the n equal to see whether this information can fas 
obtairisd from some combination of the® we therefore consider the 
linear comb inat ion 



i ffismber of 


wi thou. 



L ' 


— + 1 


( 1 . 4 > 


The condition for i 


i.2> to be in the fora* (1.4J are 


i 1 -1 


i i j j 


sliaiinate the ffi to Qive 


1 £ A 


J' 3 = 0 


< 1 . 5 > 


i 1 J 


There are h equation 


for the muitiplyers 


and the direction 


they a.r-s homopaneaus in the 1 


necessary and 


sufficient condition for a non 


trivial solution is that th 


de term! fiat ® 


coriditiors on the direction Qt 


:urve 


;o bs a. cliar 


:teristic and the corresoondinq equation 


hs in . charac tenet 1 


tic form 




Q r OP ao a 1 1 on o 


a, .sLirface exoiosion o 


f strong shock i^iaves in riait space dti© to 

r impact, has been treated in fflany iaveles 


jxunation . in one of these an attempt in made to account 


ot aDoro 


for the Biaterial strenptn by 


includiriQ a r4ew-tQniori viscosity 


fcerm. 


This approximation seems to have been oripinated bv S.H 
co-workers V291, Who have further approximated 


tf'uan and various co work? 


eouation b'- 


seek 1 no 


t hess so i u t i ons 


ordinates 


space 


of quasi "*■' siiBi Sioiut'ions' » 

.tions the Viscoitv coefficient is taken to be 
function of time but independent of space co - 
in addition , the flow field in the half 

represented as throuph it were one half of 


iPhericai iv svaimer 


flow » The predictions of thi': 


theory 


for certain quantities 


on the sym.metry axis at 


the shock 


front are then composed with experiment and it is usual Iv claimed 
that the apreement constitutes a validation of the model 
used . However it has been pointed out repeatedly that 
auantitiss a,t the shocks front are very insensitive to 
= model and can be aproximeted fairly well by 


luantitiss a,t 


t hi 


jetales of me 


any numbers of approximations . ft recant 
of this fact is presented in a paper by 


L 1 lustration 
Bill inqsley 


whose result show 


the aiffsrer 


IS that experimental data was incapable of' 
nee between a wide variety of theoriticai 


treatffients S'- 


the subset of quasi-simi lar solutions witn 
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3Ppr'D>t iiTiS't ion ?jssd m a 


definino the equation of state to mention 


ilv a Tlis 


V been cr it 


ically reviewed in a recent survey 


paper 
added t 


bv Ras 


spherical iv symmetric maqnsvit 


tic field i 


' i sccsi tv term 


Lslatiort along with a piewtonian viscositj 


Tiie above paper propos 


r.es to study the composite effect of aii 


the approH iifiat i on 


V on -Me uiTiao n an c 


Richtinyer £323 developed a ne« techni.cio.e ol 

.tv and introduced an additional term in the 


foliowino Torm 


where 


;< the internal 


dimension 


s-ionless Constant ,V the material velocity 


enoth and po the density at We have shown in 


this thesis xn 


discontinui tv in 


dinary gas ayn 


that artificial 


ss i t 'v can s mea r 


out shock 


ffiedium foore easily than in o^ 


(iii yi THAU’S RULE 


Chester C331 has studied the motion 


of a shock wave down a 


non uniform tube on the basis of a linearlized theory in which 
the changes in the tube area and the consequent changes in Bhock 


strersgth have been assumed to be smaii. in tni = i ineaf i -ed 
theory, the solutions breaks down when the flow behind the shock 
is nearly sonic Chisnel E343 and Witham £202 have shown 


wave is near! y 


sonic 


that Chesters work could be simplified and extended one minor 
simciificaticn is that where as the Chester worked with the fuii 


three dieeneiohal aquations and performed and averaginp process 
in the course of his analysis. It is sufficient to work from the 
outset with the one diamensional formulation . Witham has shown 
that the motion of the shock can be found in a simple wav without 
solvino the equations for the flow behind the shock in detail . 
This nisthad cars be givin as the faiiowino ruie . 


The relevant equations of motions are first written in 
characteristic form . Then the rule is to aopiy the differential 
rsiation which must be satisfied bv the flow quantities aiono a 
characteristic to the flow quantities just behind the shock wave . 
Tooether with the shock relations which e5<press these values in 
terms of the shock strength . the rule deter 0 iines the change in 
the Shock strength . In fact it is found to the accurate even 
when the total changes in the shock strength are not ssnall 
although a full understanding of this fact is still lacking - 
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magneto 


CHAPTER - n 

radiative shock wave propagation ih a conducting plasha 


1 . INTRODUCTION 

The problem of propaQation of shock «ave in a non 

has been studied by many authors where the 

homogieneous iisediuffl nas 

and atmosphere of the earth are 
effect of eKpiosion in the star = »na a a mu v. 

discussed . 

, A.- «-f rvlindricaiiv svfflffletric flo« has been 

The solutions of c> i mar icai a y 

■ ^ H I in r^3 Rav t23 has discussed the pr-oblems of point 

obtained by Lin 

-srtH fr-'ind an J^nact analytic solution a 
and line explosion and found an -t. 

onalvtio solution in the three oases of Plane ,cylindrically 

symmetrical . Spherically symmetrical flows have also been 

^ nw Sal-ur-ai C:^.3 . Roqer'S 1143 has also studied the 
discussed by tsaKur ai s.-.* - 

similarity solutions for ail the three cases in uniform 

atmosphere -Later on Sinph and Vishwakarma CSl have discussed the 

•isrif solutiohS of the flows behind shock waves in a 

total enerov iricrsaBes 

radiative magnatogasciynaiiicB lo which 

^ith 

In this chapter the problem of explosion along a line 

-i ft*' laritv solutions are developed 
gas cloud has been discussed , ^imAiarify 

describing the propagation of a cylindrical and sphe. ical 

in a non - uniform atmospher-s where magnetic effect i= ob=er\sd 

p„e to iOhisatiOh of pases and the problem is solved , takinp 
counter pas pressure and radiation heat fluk in to account . The 
radiation pressure and radiation enerpy are ignored due to weak 
ionisation of pase s ■ Th e pas in the undisturbed field is as su^d 

* Published in Astrophysics and Space Science (Netherlands) 




-if-, r 


2G 


to bo at root . «e also aosuoo the oas to bo ,roy aod opaauo and 

the shock to bs transparent and isothermal. The total anerp. of 

constant . The Runqe Kutta method has been 
sKplosior. IS con=.an^ 

j_ 1 at X iT kaopot" » 

used for aPBrOKlmations on computer sy.tsm Ctb 1 - . - 


<2) SELF “ SIMILAR FORMULATION 


eqiiat ions 


of conservation of mas:® . 


The fundamental 
ffiomenttsm energy and eguation 


of the 


maanatic field behind the wave are piven as 


ae 


pii 




4 - — --ipuJ J 

at Sf ' 




at 


Cpu) 


a . 

} 


s}h 


ap Sh 


< 2.21 


ar 


at 


' tpE 


) + 


3 >-' 


[puE 


) -f" 


— CPu> t 


1 J 
Cqr ) =0 


C2.3> 


dr 


ah 3 
at ar 


\^hu 


12.4) 


fe’hsre 


.j=i ,2 for cyiii 


Lndrical and spherical respect iveiy- 
and ^=i. r is the radial distance from the line of explosion, 
density p , pressure P, radial velocity u 


internal energy 


E, and h is the magnetic field. 


heat flux q, time 
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For an ideal gas internal energy E is given by 

P «=: 

_ 5 E = E (T) , 

(T-i)p 

where -r is the ratio of two specific heat . Furthermore , the 
eauation of state of the gas is taken to be of the fora . 


r p T 


i2u6i 


wliere F is tiis constant 


It is assufiied the local thersiodyaamic eoui librium and taking 
Rosseland's diffusion approximat ion CZelodvich Raizers C63>. 


3 dr 

where T is the absolute temperature !? C / 4 is the Stefan 
Boltzmann constant, C is the velocity of light and p , the mean 
free path of raidiation is a function of density and absolute 

tefaperatiirie I « 

Foiio^^inig WBnq £7i 


« P 

—• fis ^ T ' ' 

0 

where |i , « and fJ are constants .The disturbance is headed by 

O 

an isothemai shock surface , hence the conditi.on across it are 
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C3) SOLUTIONS OF THE EQUATIONS 


In front of shock 
density , pressure and 


in the undisturbed passous mediuffl 
maonetic field distribution are 


7 


tiie 


p = AP! 

i 


—2 < M < 0 


C3« 1 1 


n 

P = BR 
1 


C3-2I 


II 


ii = DF< 

i «i ® 




where R is the shcok radius at the tiras t 
H, are constant . 


A , y , B , n 


D and 


!he 

governing 


similarity 
the flow to 


variables which reduces the equation 
ordinary differential equation is taken a.s 


n = 


a b 
r t 


f3«4) 


where a and b are constant . It is assumed that the flow is 
headed by a shock front given by 

0 = 0 , 

c 
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the other tranefrometion for the flo« variables are 


r 

y _ VCilJ 

t 


t3«5i 


k X 

p = p tiltisJ <! 


hf 


k+: 


p ifO 


(3.7) 


(k+2) ^ 


h = 



t 


H (ii) e 


(3.8) 





F cri) 


C3.9) 


a b 

where O = r % e 


The total energy of the disturbance perunit length is 


R i 2 

q = 2 kJ ft — p u 


P h 

^ + 3rdr 

( T - i ) 2 


(3. iO) 


y 


6 


Intersus of %'ariable O we have 


*"^0 <-> p 

rs-yr I ? ' ' “1 ^ W' 

Q = ----- J i: - ^ 7 — :: ■" ^ 

a O 2 


a 

<k+4j / ia -1) C A- 2 — <k+4> 

«■ n ^ ^ 

LOhcTt Hq Hie Kjtxkxx^ cstr n °-t 

n = E , o *N -^ = 

where .Shccksurf acs to be given by n ^ = 


* dn , <.i, ii) 

Hie v^h<rck i-rsni 

■■ 1 

constant . This fixes 


the 


velocity of the shcok surface as 


a t 

^ #"r^' 4 r*irk <KEii''''f jsr* ® if a i • The tot a 1 en e? gy 

which represents an out going suriace it a 

of disturbance within the shock surface at ' any time tis 
coristant » This by C3»ii) requires that 


S 

^ _ 2 — - Ck+4) = 0 . 

b 



Let the Mach and Aifven Mach number's at the shock front 


l^^gpect ivsiy ias dafioed by 



ion of (3.5) - C3.9> in ths equ 
2 . 1 '} - <2.S) shock condition (2.9) - (2.13) and after 

■elation (3.13) >j find the similarity condition ^ 

inqh and ths Vishwakarma £53 are cofiipatable when 


utilising 
f oi lowing 




w + 2 

and It = ^ ^ 

tell c e p SLi r*© d i t r i Is ii 1 1 on Is sc omss 



C3-1-?; 


P = BR 
I 


(3. 16) 


and equation (ii. 1) 


(2.4) and (2.7) are now transformed into the 


forms 



34 


H'' cn) 


H CO) 


1 

V' <ri)= — 
H 


P CPi) 


fHO) 


The appro} 


V<ri i 

0 


iicn > = 
o 


ni4+m v'dl) - tSsi/2 + 0 + 2 ) V (ID 
n i;2 - <4+i4) V '(fi) 3 


h. 


T v*’7' ^ 


a 

(F ifii/Pii i fi (ri)/A> C' 


V (11) 3 




i /2 


L P(fl)/A 3 


r!tT2 — (4+wi V tfs) 'i^CA-f"^) 


<J+1) C4-^M> V in) 3 


E2 W + 2 (t— i ) H <n) 

n C £2 - (4+w> V Ul) J-C44-UJ) j 


2 . ..o^riA,«,_Ln EilLl 


t2 ~ <4+W) 


Vin)3 Vini •£ V in) <i -J) - 1 3 


(3.23) 


n E i: 


(4+W) V (o'* 




lissate shock condition are 


£ i 

C4-t-w) 


1 i 

2 2 
T H T fi 

A 


(3.25) 


7 II M 


A 


(H + H ) 
A 


iO 


(3.26) 
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P <0 ) 
0 


4 


C 4 - 5 -^> 


H 


A 


( Si + M > 

A 


(3.27) 


j“ \ I I ) 

O 


2 


H + M 
A 


^ 4 

(4+y> C T !1 + M 3 

A 


(3. 28) 


H (n ) = 

o 


4 


2 T H 


4 


A 


C 4+1*1) (T M ) 


H < H + M ) 


(3.29) 


i^ihen ll — i 
0 


C4) HESOLTS AMD DISCUSSiOM 


For BKihlbitinia the numerical solution it is convenient 

the flow and field variables in the following non - dimensional 

forms 


y 


1/ (4+w> 

<n ) 
o 

(O) 


y (O) 

V (fs ) 
0 


C4. 1) 


m/ < 4 +w> 




Si] 


?_ 

p 


(O i 
0 

tni 


Q (ii) 


si ifi > 


< 4 » 2 ) 


P 


Cw+Z) / <w+ 4 ) 
if! ) 

o 

cn) 


P (I!) 

P (0 ) 

o 


< 4 . 3 ;* 


(w^-3) / Cw+4 ) 




cn ) 

0 

ill) 


F (fl) 

F til J 
O 


<4.4> 


C I / 2 C ) 


h 

h 


<n ) 
o 

cri) 


H <n) 

H tO ) 
O 


t4.5> 


Numerical results for certain choice of parameter are repoduceP 
in a graphical forms. Calculations are made for folloeinp value 

of parameters. 


C i ) ^= 4/3 


=20 e M =30 , W= -1-5 


ot*i/3 , 


j=i , 14=10 '0=1 


A 


< i i > ‘t - 


=4/3 , M =20 , n =30 , W = -1-5 , «=!/-:> , 


J=i 


0 = 1 , J=i ? 


A 
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=ioo 
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C ii i 


> t= 4/3 , li 


2 2 , ^ /T T~<? 0=1 ,N=iO 

=20 , M =30 , S^= -l.-i f ^ ’ 

A 


^=4/3 , H =20 , =30 , 


W= -1.5 I, «=l/3 , 


j=2 , ^=1 . N=i‘>--’ 


. J4. cls=>ar that the velocitVf density 

From the sjraphs it i= cleat 

pr«s.ore,™a,n.tic fiald distribution and heat flu., are »,a.,ieue 
the Shock front and decrease eore rapidly towards the line oT 
explosion as the value of fl increases and there is aiso^ Blight 
charge in these variables »ith the change of radiation parmeter . 
Throuoh the oraphs it is also clear that decrease in velocity , 
ppcaity , pressure , radiation fluk and eapnetic field variables 
is more in spherical case than decrease in cylmdr ical 
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CHAPTER in 


SELF SIMiLAR HEGNETO GAS DYNAMIC CYLINDRICAL SHOCK WA¥ES 
( 1 ) I HTeODOCT I ON 


ft cylindricsl i^iave of s^plosion with a shock surface as 
wavs frcnt oroducsd on accaant of a sudden rsisass cf a finite 
sitount ot ansrgy, 3 Kpandinp oatwords in a conducting gas susjected 
to a .'.magnetic field has been studied by Psi C i 3 ,Chr ister and 
Hslliwale £21 and sany others. Dsb Ray £33 ano Singh - VishwaKarse 
E43ha'/s oir-tainsd siiTiiarity solution for strong cyiinoricai blast 
|s; 5 ave'S in a conduct i nig non Liriitorni iiediy.iTi *» 


Th-S present ctiepter deals tfis pnobisin 


of the 

instate 

neous 

re 1 ease 0 f en s rq y 

a !■ OT"i eg a. i x ri s 

in 

a §'§.s 

•C 1 OLld 

thunder 

in tiie effect of 

mapostic: fisid 

where 

density 


with 

time*. The motion of b 

pas is as'SCifiisd 

to be 

adsabat 

ic and 

tota, i 

ca:|'"5Cft6-'’-Q V of the i*j3,ve 

rsrsains constant. A 


cofiparison has been fuads on 
axial coffiPQnsnts cf the 
sunfsicS'* The tr-5.nsvsnss and 
’'1“' ? I cs !*: 1*^' f •^*ns h sod h nesp 


the effects of the t nans verse and 
ic field varialiies Dehind the shock 
axial cO’iBpanents of the Tiegnstic 
actively » The oLisierical calcuiatscn 


ha 


on Dt 


system 


R. 


, 6 s S d roar am® = 


C2J EQllATiONS OF WOTIOli hMQ BOUMDAKY CONDITIONS 


Fcilo. inc yithsifi LdJ « trie eouatiGn oi iiics.y i 1 5*5^5 

«n£s-™:ri "i "| 3 ' 3 rrr i 1 ansteadv of a perfect pas irixtn transverse ajid 
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^ H 
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a” 


1*5 
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/ T "T \ 

^ -mJ' • Is-#* i 


Ths dot 


t i 'iiiS arid p 

rifiie 1 

sifs5i lar-ity 

var i ^ 

eaL'.ation be 

comas 


:istion of scales with rssosct to 
;itian of functions with respect to 
isinp the relation, dif fsrsntial 
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f? 




C ¥ •" III 


3 = O 


C3. 4) 
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•i -I- Cv - HI 


2 p 
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C 3 . 5 I 


R R 


L V 
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1 1 ' 


) 3 


C3. 6> 


o 



Inorder th- 


iniiaritv solutions C3»kf)' fee near 


necessary that variable t and ii in equation c; 3«4 


rstist ba 


separable 


:his pQssiois It 


Lind p 


ions •k-sn > 


jra ft , B, « and 8 are constar 


shus aii the 


in the 


self similar mation have a cower law deosrd 


srdsncs on time . By use 


relation (3.9) and (3.10) eauations (3.4) 


<3.S) take the 


non d ifnension.: 


srfp of ordinary aiTTs? 



t is clear fro® the equation 
wavs takes place »liari B = “- 2 « 
siqrii'^ic.ance only when - 2 IS 

A 1 1 s r os i fi q s i m i i an t y 


C3wl9> that unifor® expafision .of 
The BGiution have thsir chysical 

< o 

transfomiat iori C3*2) tiis jtiiiip 


’c /) trscooss 



h 

e 


whs re II 


(4) RESULT AMI) inSCUSSim 
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?+•! 


"f Ft I 
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' 4 ^ 4 > 
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T-^l 


H] 
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T --1 


f1 h <11) 
ft 


zi 


Numerical result for certain choices of parameters are reoroduced 
in the form of graDhical form . The caicuiation of the result are 
for follawinG values of oarameters . 


< i ) T= 


s= 0-25 , H =10 , 
A 


=i. 4 


«= 


N =10 . B=0 
A 


in this arcoiem. for the shocv; ^;hic 
eKCicsicr fncves ;.r cari'ci.ictr.no cas-we have 
mayes a cone iri aicr. csTUSsn than . ' 


h is oroaucad bv a line 
ci f ctissed tKC cases and 


Css 



It is 
a = 
and 

surf ac 


cc 3 ®=arisiQn has also been made for both the above cases . 
observed that when the shock radius is contrating that is 
25» ther is sharp increase in velocity , density, oressure 
;=snstic fsild of the fluid particles behind the shock 
■s as compared to s( = -50 , when the contraction into 
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CHAPTER IV 



PROPAGATION OF SPHERICALLY SYMHETRICAL DISCONTINUITIES WITH 


INCREASING ENERGY IN GENERALIZED ROCHE MODEL 


(1) INTRODUCTION 


isctissed methods for obtainins^ analytical solution of the same 


non 


:Oiution in the Ben 


Roche model I 


arround a nucleus having a iar 



(2) EQUATION OF MOTION AND BOUNDARY CONDITIONS 

The equataan af continuity aioinentuffli . field and energy in the 
Bsneraiized Roche model are . 




a ^ 


au 

S r 


pu 


=o 


■{ 2 . i ) 


— + 

at 


att 


— +■ 


1 aP 

p a r 


__a 

a*" 


2 1 h 6f'! 

th > + "S' 

p r 


O , (2.2) 


ah ah 

dt a r 


a u 


+ h — 


=0 


(2.3) 


ap a p 

+ u — 

at dr 


tP 


ap 

at 


ap 

+ u 


^ a 

+ ('T-l? (r.F) 

r ar 


U 


(2. 4) 


«ners V- ,P,^.h and F are the velocity , pressure , density , 
aiaqnstic field and radiation heat flux , at radial distance <r) 
froffi the centre of core at time (t? ? 6 be the gravitatianai 

constant .The magnetic permeability of the medium has taken to 
be unity through out the problems . The equation of state for 
ideal gas is given by 


P = rp T 


(2.5) 


where F is the gas 


constant 


61 


Also assuffiioQ 


local thermodynaroic equilibrium and takinq 


Rosssland's diffusion approximation iZeldovich RaiEors t6> 3 


Cji a 


hr 


4 

<cr T > 


( 2 . 6 ) 


where cfC/ 4 is the Stefan - Boltzmann constant , C the velocity of 
lipht and |i the mean free path of radiation is a fanctxon of 
density and absolute temperature T . 

Followinp yanp £73 


oc B 

= p. p T 


(2.7) 


o 


where p 


. fin beinq constant . In the seif similar model the 


total enerqy of the wave is dependent of time as 


d 


E = B t 


(q < 0) 


( 2 . 8 ) 


where B and q are constant 

The flow varisblBO i»,medi.tolv .head of ehock denotsd by euffXK 


( 1 > are 


“W 


U = 0 , p = AR 
t i 


(0<W<2) 


(2-9) 


whoro ft id a constant and R danotB. the radius of the ehook 
eurtaoe ahead the shook , the oagnetid field distribution 


IS 


h = CR 


2i = W+1 
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and pressure distribution ahead the shock . 


C 1 -J-w I 


A6HR 


<1”0> 


p ^ 

1 imii 


2 -20 
C R 


(2- ii > 


where C. ^ and 0 are 
conditions headed bv an isothermal shock is 


constant . The Rankine Huoaniot shock 
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C2- 14) 
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u > = h 
2 1 


(2. 15) 


f.2» 16) 


where suffix 2 , denotes 


and 1 denotes flow van 


Ifl as 5 perLtn 


it area across the 


the flow variable just behind the shock 
iable just ahead the shock denotes the 

the shock and V be the shock velocity and 
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qiven fejy 


dR 

dt 


C2. i7> 


TBAt^SFORMATlON OF EQUATIONS OF MOTION 


se equation of flow to ordinarv differential 


in order to reduce the 

i ^ -t-hpi folio 5 .<iinQ transf orsiiat ions 
equation we nas^^ introduce xne 

Follow! HQ Sedov £73 


-- 1/3 - 0 - 

n = ( 0 £ HB > 


< 3 . 1 ) 


where ^ 


2 +q 

5 -W 


<2 - 


( 3 . 2 ) 


and the limit of d and w are 
0 <.a < 4/3 and 0 <W <2 


( 3 . 3 ) 


‘ 4 a sea 


the solution of equation 


< 2 . 1 ) - ( 2 . 4 ) in the form 


— yen) n ^ 

t 


AMS t 

W+J 


Rcn) 


p 


AHS 

W+1 


ptn> 


i AMG > 
w+i 


hCfl) , C3.4) 


AHe 

w+i 

r t 


FCfn 


Jisincs eaustion t3.i) in the equation (2.&) v^iith the help 
of eauation <2. 5) we obtian 


OL 


W 


w+ i 


and !J 


-(514+7) 

2(W+1) 


The equation (2.1) - (2.4) and the equation (2.5) are then 

transformed with the help of the relations <3. i) and (3.3) to 
foi lowing farm 


vcn) 


f3~a+4 __ ^ 

1 FCfDP ^ , 

( V-cr) (V-ff ) i:2f lp+ (i4+i ) (p+H .* -H K 

.1 ’ ■ <55* * 

N t*+--* 

Ptfi) 


R (1-W) 

1 + 1 44+ 1 ) V*p- 2p VH 


LfiP -flR CV - ^ 


where 


4 Ciiosr ® ^ 

— (AMS) 

3 4+0 


(3.6) 


is s 


dimension less parameter 


6 
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f 

sen) 

RCO) 


1 

CV-0^) 


(W+i) ' 

VCR) - V (0) 

n 


n 


(.3. ?> 


/ Cw+l) 2 H R R VRCV-t) 

p (O) = EP + H 3 - 

n 0 4 Ff 

a n 


H H , <3-81 


i i FlRP > 2tP 

Ftil) = C P 3 + 

(t- 1 > R 


‘Lw < T-'l ) +3 t - 1 3 


<T-i)n n(T-i) 


i 1-k) 

n 


F <n> 


<3. 9) 


4. RESULT AND DISCUSSIONS . 

The transformed jump condition at the shock front is given by 

2 i 1 

v<i> C i- 3 , (4.1) 

( 4+1*# ) 2 2 

Til! TJ13 

A 

*7 7 

tH m 

A 

Ri;i) = — , (4.2) 

2 2 
< M + H ) 

A 


7 
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p 

p 


1 

14+ 1 


P(f» 

P(i> 


C4.8) 


n 


h 

h 


i Him 

iifi+l}/2 H (1) 

n 


<4. 9) 


1 F tnj 

« F<i> 


{ 5 . 0 > 


n 


The numerical integration is 


carried out on DEC-system 1090 


comouter installed at I.l.T. Kanour by well known 
proaramme for the values of w=i,l-5 . 

The other parameters are 


R. K- S, S 


n 


15, M -25, a=i/3,S= 


• 1 7 +5« ) 


2 < i+w> 


H==io 


The nature of flow and field variables for adiabatic case are 

illustrated throuoh graphs, it is clear froe the graphs that 

s * X. HI I n i fflufii 

velocity , density and pressure 

shock front but increases rapldely tpwords the center of ekplosion 

in adiabatic case . Sphere as distribution of mapnetic field and 

-+• ^hoct front but decreases 
radiator, heat flux is maximum «t =hocK 

towords the centre of explossion. 
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CHAPTER 


*NALYS.S of self SIM1U.R H0T.0H in the th^ 


EXPLOSION 


(i) iNTr'RODUCTIOH 

car, -US stal 111 .nO Sedov C23 were first to 
O.scuss toe eodel of etell.r etoiosioo in whioh a star is 
considered to be . perfect self Qravitatinp oas. The distribution 
pas at any eooent of ties spherically sveeetric. In the 
,_Pnt Chapter a anyalytlcal solution of the classical eodel of 

atellar e.plosion has been invastipated. . nueber of new solutions 

1 1 j»tion of Qas occur after 
has been obtained in i^hicb radial oscillation 

. Takino Newtonian gravitation into account 

the sli Q'C k wa. ve p as »es « i a k 1 1 1 u 

a. ^ ^-e -J=lf similar motion of a Qss 

a thorough analytical s y = 

„Her the effect of mapnetlc field developed in the 

dynamics under the exTsui. ... 

■ rs irfF-irh was earlier applied by 

theory of stellar explosioh which 

,=ny..i:4l It is supposed that originally the 

NQvikovC33 and BogoyavienkyE:43. It 

1 -he oas density p , the pressure P, 
star is in equilibrium state the gas d 

¥Hin sphere of radious r ^the radial 

the mass M of the gas within P 

a. 1 .. field h have the form 

gas velocity lu the magnetic xieia 


tla i> 


Cr ,2k*w+i. 
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2«A 6 


(3-w,f Cw-1> 


2Ci~w) 


C (l-k> 


where A.w and 6 are constant. As a result of enerov 
libration at the centre of svmmetrv r=0! a shock isave travels 
fPQS, the centre. The motion of oas behind the shock Tront seif 

similar and adiabatic 




78 


Sh 5^ 

Bu u 

4- u ■+ 

h + h — 

Bt 

9r r 


2 

s= 

4wr • 

d r 



o 


<2.4> 


t2.5> 


^ere are radial distance from centre . time . 

veiocitv .densitv . oressure and mass contained in a sonei e of 

rad i us r 

The disturbance is headed bv an is 


othermal shock with condition 


P. 


<v-u ) = V = m 

2 1 s 


(2-&> 


h 


p ^ — o 




m u 


i 2 


< 2 . 7 ) 


h (v—u J == h V 


< 2 . 8 ) 


t + — 

2 p_ 


1 2 2 
— (v-u ) * 

2 2 


1 1 2 

c 4 - + V * 

1 ip 2 ^ 


C2. 9> 


C2.10) 


T = T 

1 2 


C2. 113» 


n = M 


„here suffix 1 .nd 2 denotes the flow variables iust ahead and 
"ust behind the shoot- front respectively denotes the mass 
flux oerunit area across shock and V denote velocity 

Shock , is ginen by 


dr 


dt 


(2. 12 j 


a , tr Sedov C 5 1 . the total eneroy of pa^ between 

AccordiriQ to qsgov l «? 


spheres of radii r , and in eaui I ibriu® is 

i 2 


5*^2w 


. 1 ^ 

E := 3i^ SA 


1-2 {w-l> r 

(t- 1> (w-i5 (3-w> (5-2wJ 


(2. 13) 


3271 {4-3t) 




1 oq 

r 

1 


cr tartf - 1 r»t 5 *=»d wiiihin ■fehs spHet s 
If w < 5/2 , the total eneray E enclosed w 

, ,f „ i 5/2 it iS infinite but E < u ir 

ai radii r le tinite i r « - -i- 

r*i I /’ m liie 

, > and E 1 o if f i r where f =--i ^ 

o 

. 1 -s , Eiv'nl nsions for w 2 5/ 2 in the 

- *■ - the nrobiem of stellar explosions 

solutions of tne pTuuif=« ^ 

1 sDoliscI to 

frtrtfion considered can only be appu^f 
class of qas motion cu«ta 

e - „.vr.=i-.i-e of E at the lower 
, hiacause of divergence oi t - - 

astrophvsical solutions because 

1 are regarded as intermediat asymptotic 

limit . These solutions are r ega 

out Side a small neighbourhood of r-o, 
solutions valid QUt. ^ f=+.-hio 

. , (E . o) discribe the break down Of unstable 

sciiitioos with T -i ^ 

o 

.teiiar eouiiibriuo ■ - iaw of enerpv iiberatinn in a Seif 

,t.iiar solutions has the fore i Bopoyasienky t4i > 


4 
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<5-w> /w 


2 115-2^3 / w 


Evidently E is independent of the time for 


C2. 14> 


5/2 or ffi==0 


The constant a is calculated frora the solution it self and in so®e 
case turn out to be infinite . The corresoondino solutions then 
provided asvfiiototic form for a very intense explosion . The 


faliowioQ ne«»i resul 


(1) For T 


ts mav 


directly be deduced bv 


. = 5/2 and also for t < 2 iw-l> / 3 there 


.Kist no solution with a vacuum formino withen the oas 


For 4/3 • w “ 5/ -a- 

accur after the shock wave 
limit of dynamic svstem 


M=i . damped Oscillation of the oas 
passes which are connected with the 


(2) For T < 4/3 , « “ S/2 *11 solutions have a soherically 

vacua®, of increasino radl. which fores aoout the centre . the 

oas iHonotonicallv spread i no out from it . 


[3> For T < 1/ 


- the oas returns to eaui librium 


after the shock wave passes but in case when 


3 + C2w-5> 3 


3(w-i j 


repeated damped oscillation taken place 


SiHILARlTY SOLUTIONS 


The Similarity variables which reduce the equatio, 

flow to ordinary differential equation is taken as 


equation govern inq the 



Rim 


6fc 


7f) 


p = 


PtO) 


4 

t6t) 


H = 


nifo 


6t 




( 3 - i > 


- l/w 


where ri = r (AGt) 


Bv us 
as 


.„a relation (3.1> the differential eauation are 


transf ormed 


2 3 ‘’^ 

V + V + o ^ 

dll 


2 2 d 

3 - + ■ — ~ — + 

U r R 


n 9P 2^^ 
in R 


N 

+ H = O 

R an 


< 3 . 2 > 


an 

+4p3 V m>!;2p + n --- 

M dn 


TP 

3 t-2 I - 

R 


n 3 R 


BR 


4- R 1 +0 V 

w an 

1 o e C3.3> 
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..pstioate nuiiisricai flQ« solution we wrix^e txuv. 


Inves 


diaiention fom 


y <05 


ROD 


R < 1 5 


p<ro 


p<i) 


h N<n) 

h Hil> 


M 3 fieri) 

M * M < 1 ) 


ISCOSSIOM 


Above j-'^eiattons 


the dietribution of velocity 


„ maonetio tibld and maes distribution in 

.sity, pressure , ^ ^ 

sllar model , when magnetic tie 

iliar moae illustrate 

, also be obtained which may 
3ro'.'ii-nation may also oe 

. . pvactlv beheind the sur 

n -,an,r! field variables exacci, 

laviour of 

discont inuity. 








CHAPTER - VI 


■HOCK WAVE OF VARIABLE STeE^4GTH 


Ml EXACT SOLUTION OF NORHAL 


ADVANCING INTO A REGION OF VARIABLE DENSITY 


I NTRQDUCT 1 ON 


workers to 


accord ino to some inverse sower of the 


FUNDAHENTAL EQUATIONS 







where 
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f^' 






U = 


< 2 . 7 > 


3t 


hP X 

and also = a (ps+R^) 

as 


< 2.81 


Now,wa try to find an e 
satiations of leotion^ 


Kact and special solution of the 


( 2 ) SOLUTIOM TO THE PROBLEM 
We choose a variable 


{“ks/s^ 


< 3 . 15 


n =t a 


where « is a 


constant to be determined. To obtain the 


solution, we can take the van 


abies in the foilowinp form 


1 ^ 1^ cn 

;,;=t xf.n5,u=t u<n),e=t iicn> , 


C3.2) 


2.ti+c^2. 

p ' Pill) « 


:t to n , equations < 2 . 55 , < 2.85 can 


Biffsrentiating< 3 , 2 ) with respec. 

^dinary differential equations for i,U,i 2 ,and p as 


be rsdiiced to or^ 


I i ^ 

: CC si X + fi O X + n 


Q U5 li si K = d , 


( 3.35 


/ f 

p + fi X lid 


I ) u + n u 1 ~ o 


( 3.45 















where A 





reduce to following form 


Y + ic+a)Z + Y <b-i> - Y 2 


[3. 17 > 


where H = 


ih-i ) <T-1 j 


— Z C2C+2b+a) +2b+L;— T (C+«) + 


+ <1-T- 


which atlhough not very simple but may 

In order to find the quantities 

_i fnt- V 8f 2 from <3.17) 

<3. 1?) by quadrature lOf i «£ ^ 


determine the solution. 


WB have to solve 


we derive 


1 dfi 

n dz 


;z-i> 2 


c (T-.1) Y 3 + (2b+2C+«) 2 


(T-l)+T.a-2b 


iiT-D+Z I > ■ 


(3. 18) 


The shock problem suggest 


s that c 8. « is such that c+«-0. Thi! 


leads to a simp 


iification of 


(3. 17) . The va 


lue b=i corresponds to 


tti® cass o 


,f homentropic flow. 


4.STS0I4G SHOCK COHDITiOHS^ 


- 4 -H a -eoion of gas which is at rest 

It the shock is eovinq «tth a .eQtc . , . f , 

„ but a un-uniform temperature t. density. 

under uniform pressure.pf 

suppose that density distributinn obeya tbe pdeer ioe 

, .-..<'"„her. t t. e are constants and the suffik 1 denotes 

e. p,Mtion at shock at the 

the Shock condition ahead the shock. P 









Putting 


fr'‘OffiC2.4) toC 
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CHAPTER V I I 



analytical solutions of cylindrical shock waves in a rotating CAS 


WITH AZIHUTHAL MAGNETIC FIELD 


INTRODUCTION 


Kumar and 


and diverging 


through a rotating ga' 


imuthal 


i mu i t aneous 1 y 


i distribution is 


and shock 


hoci: 


the magnetic field is weak and when the 
■or strong shock also we have considered 
inetic field is strong and when 


ofidit ions 


non 


wliicfi xs punsi/ 


case* 



(2 > BASIC EQUATIONS BOUNDARY CONDI 1 


lUN ANU a™ 


for shock velocity 


the cylindrical! 
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muthal and 


ihack conditions can be 


ramediately 


hock 



9.9 




WEAK SHOCKS 


For a very weak shock we take the parameter M as 


e 


N=1'*'€ ,E <<1 


< 3 > 


y 


Now consider the two cases of weak and stronp macnetic field 


CASE *• i 


Whan 


the fflagoetic field is weak i*6« b 


o 


a iimBm |iH 'vtF* ) 

o ■ 6 0 


and therefore 


the boundary conditions C2> reduce to 


p = f t; 1+ €) , u 


b 6, H - H 


< 1 + €) 


y 


©1 BQ 


H, =H, Ci+e),P =p 
i 0 i O 


[fe) 


T + 1 

|j=ti+ ^ 

4 0 
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CASE II 




When magnetic field is strong »a ^^(i.e 

under these circumstances the boundary condition (2> 


p i 1 + € ) , u 


1 O 


=b ,H =H (1 + fc ) 

i o 01 00 


H =H 
z 1 zO 


a + E u p =p ( 1 + T e ) 


1 0 


U=!(l+ — 6) b 
4 0 


STRONG shock: 


CASE - I 


When the magnetic field is strong 


b >>a 9 

o o 


under this 


condition the boundary conditions 


a 


) 


reduce to 


(7> 


2 

pH i 

o 0 


(2) becomes 
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2 

q tz u 

dP+a<H dH +H dH > + pcdu + C 

$ O s s u+c 


|[iH (u— c) ' 

0 ^cv 


u+c 


U+£ 


dr 

r 


0 , ( 10 ) 


.■” 7 * 

where C = a + b 


rp 


H + H ) 
6 z 

p P 


and H = H + H 

o # 


WEAK SHOCK WITH WEAK MAGNET iC FIELD 

i9) and nectlecting the 

Sufast i tut inig shock condition 

second and higher order term oT fc , =.^nc_ 


dP 


da 


o 


( 


H + H ) } dG + '1 


ZO 

O 


P 


2pM 2dP 

rir ^ 60 

tP 


Q 


o 

K 0 


} e 


tP 


(11) 


dP 


da 




Q 


fiH 


6o 


:o 


Now substituting 


and 


tP 


tP 


O 


a 


o 


o 
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in equation ill} we get 


de 

€ 


ET(i-4e 


T0 


IJ 


A 


dr 


log 


( 12'J 


■? 


where <3 


U<H + H ) 
6o 2Q 

T K 


11 H 




T K 


Integration of equation (i2> yields ? 


€ =K r 


tTti-AU >-43/'r,e ^ 

'T‘ 1 eKp t ■" — T J i ? 

20 ^ 


(i3a) 


where K is constant of integration . 


In the case of diverging 


weak shock with weak magnetic field, ws 


have the relation 


JU4 


-1/2 ft -t) 

e =}•< r (loQ ) i exp' 


3 2 A ~i 

$ (log ) 

8 i r 


€13b) 


where K is c:i3nstant of integration and 


q =(■ 


hence from equa 


stion (13a) and (13i3) we have 


[7(1-40 )-4)/T0 
- 1 1 expf.- 


•(Iqq ) 3 


A -q 

Clog ) 1 exp 


3 2 A “1 

0 (log ) 3< 

8 i r 


3 r converging and diverging shock respectively 


yEAK SHOCK yiTH STROHG HAGNETIC FIELD 


j.bstituting the shock conditions 


(7) in (9) we get 






in ecfiiatiQii <i6) 


where 


JOB 


K is constant of integration, eouation <18a> and <18b) 


4 


disclose that 


u 


A -q 


— ji; <log — i 


( 19 ) 


4 4 


O 


A -1 


and ^ 

u 3 ft e -i 2 2 

= 1 + — K Clog — ) ^ exp£-2f} (i-{3 )(1qq - ) ) , 

4 4 r i 1 


( 20 ) 


bo 


for converging and diverging shocks respectively 
STRONG SHOCK 


Substiuting shock condition (8) into (9), we get 


dp 

2 X N-i 2 X o 

dU £ — 3-i-U I 

Y 2 N 


da 


Y P 


o 


o ;<(N-1) 

dr 


r 

o 



dp 


2 2 

tk; n 0 

a i 1 

f P 

Q Q 


<N-1) + 
CN-1) - 


4 N!r 


dr N 

r C(N-l) 


3 

( 21 ) 


3»0, 


Substituting: the respective values , we get 
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u 


1 _ _B C 

2 2^5 ® 

a K Log 

Q 2 r 


C 25 ^ 


U 


A 


-B C 

< K r + > 

5 B 


(26) 


K Log 


respectively . 

■4-i^£s wsa 1 nr" i tv « tHe geosity^? 
in trie last, the sKpressions for tns velocit:,/,-^ 

^ t-t iiict behind the shock surface are. 

and the pressure of the 9 a= 3U-t 


(1 ) 


Ueak shock with weak Hagnetic Field 


^ Lr (i-40 >-4>/rff 

u = k k (log — > •" ^ ^ 

« 3 r 


I A 2 

. — ( loQ — "^1 C'27a.) 

2 ^ 

20 

i 


1.1 


k k t"" 

•rpi, 


A C i, " 

(log — > 
r 


iq } 

1 eKpf- 


2 ft 

0 ( log — 

1 r 


{27b) 


E7(i-4.a >-43/T.e 




£ i-t-k r 


i e5<p £- 


1 A 2 

— (log — ) ^ 


2}3 


(28a) 


e “ 


/ A -q 

Cl+k r (log — ^ ^ eKP t 


3 2 A ""1 

— 0 (log — ) 


8 1 


(28b) 





hock With Strong Hagnetic 


IW 



(3) Strong Shock With Strong Magnetic Field 
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CHAPTER Vn I 

effect of AfiTIFiCSAL VISCOSITY ON THE EXPANSION OF 
DISCONTINUITIES IN A ROTATING INTERPLANETRY MEDIUM 


(1 ) INTRODUCTION 

Ths propagation of strong shock waves in space, due to a 
surface e;-;plosion or impact has been treated in many levels of 
approximation. In one of these, an attempt is made to account for 
the material strength by including Newtonian viscosity term . 
S. W- Yuan£ 1 , 2, 33 approx imated the solutions of equations fay seeking 
quasi- similar solution . In all of these solution the viscosity 
coefficient is taken to be at most a function of time but 
independent of space coordiantes. Pai £4,53 and Kumar £63 
discussed the propagation of hydromaqnet ic cylindrical shock 
through a self — gravitating gas showing its vsiocity only for 
strong shocks . Recently Singh and iiishra £73 obtained analytical 
relations for shock velocity and shock strength and the 
expressions for the pressure, the density and the particle 
velocity ifnffiediatly behind the shock assuming the fact that the 
initial density and azimuthal magnetic field distributions 
variables. 

In the present chapter the characteristic method (Chester £83 
Uithaffi £93) is applied to abtain expressions of the density , 
the pressure , the particle velocity just behind the shock 
propagating in a rotating atmosphere - The effect of coriolis 
force is taken into account . Since the velocity effect has 
tendency to ' smoothen out such discontinuities, the artfical 
viscosity coefficient suggested by Ritchmyer Von Neumann £103 

■■ 'i-i-; ■ y:: 



(i^' far weak shocks and <ii) for- stronQ shacks respectively 


2) BASIC EQUATION BOUNDAKY COND I T i ONS AND ANALYTICAL EIPiESSION 


FOK SHOCK WELOCITY 


The equations governing the cy 1 indricai iy 


mfnetr 


own 


ielcl» Fol lowiriQ yithani E7J ar 


lU 


2 2 S u 


— or K 




1 - t 


where u <r,t), p<r,t), H(r,t) vrCr,t), and p(r,t) respectively 

represent the velocity , pressure , magnetic field at distance 


r)and time (t> and the density where a- 
.s the artificial viscosity. 


Tp/p 


and ' q‘ 


'he magnetic hydrodyanami 


mic condition can be written in terms 


of a single parameter N = p / p as 

1 0 


O = Np , H = HH ? u ~ il 
o 1 Q 1 


— 1 U 


y : C a + 

(T+i) - i't—ii o 


t) W +t ) 3 


2 p 

0 2 

1 Q (T-1-1)" o 


+ T-i b 


<N-ll 3 


where o 


and i respectively stand for the states just ahead arsd 


just behind the shock front ; U is the shack^ velocity, a» is the 

iw / « 'iVh and b is the ftlfven speed 

sound speed given by v , a o ^ 


given by 


% , where fi is 


magnetic permeability 




) yEAK SHOCK 


bounda 



) STRONG SHOCK 


be broiight about in 


two w 


when 


when pH 


Ii4hen b 


with 


now 


+ Yim 


H (T+l> 


Yim 




X<N) 


J17 


For divaraing shock the characteristic form of system of equation 
»■•!}, that 2s the forn in which equation contains derivatives only 


. b J p lane is 


u Or 


dP+ ^cdu+HdH 


<u+c> 


<u+c) 


CQ r 


<u+c: 


it-Ddr qtT-Du 


(u+c } 


<u+c) 


where c 


a + b 


No„ substitute the shock condition 4 or S or 6 into 

relation (7) -A first order differential eduation in (rl or u is 
obtained ehich gives the shock . For »eak shock in the presence of 


transverse isagn&tic fiexd 


on sufostifcLiting tiie, shack 


nsQierted the sacond and lie i qher. 'tems 

conditiciri 14) in 


of .£ C€ <i) 5 we obtain 
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2 

2 


dP 


dH 

H 


- 2 Q 

1 

o 

Q 

dr 

■a 

“f”- — « 




o o 

2 

e 

P 

p 

o 


Q ■ 

„ Q 



2 dr 


V 


dr - 3 . , C 8 > 


T 


Q 


But by hydrostatic eqfuiiibriuiii prevalent in front of. ■ the stiock 


2 i1 r 


dP dH H 

o i o ' o 


Q 


o 


dr 

r 


•dr. 




Q 


C 9 ) 


Hence equtian C8i .reduce to 


d€ , i 
€ 2 


0 

X 

dp 

Q 

dH 

a 

dr 



i — 

-I- 

- 1 - — 

2 tp 

p 

TP 

r 

o 

o 

o 



^ da Cl dr 

Q r q Ct— i I dr '. 

a TP TP r 

o . o 

(i 0 > 


Also we have 


-w 

pr .V 

p = - C- 

O ^ 


iw- 1 ) H :■ 

— — c; r' 3 


£ 1 ^ 


Incerting w =w/2 , ws get 


p = C 



2 

V 

H 

c 

c 

m ■ 

2 w 


w 


-3 r 


o 



where ^ 


w c w-"l ) 

Cr dq-r q(T~l)dr3 


8 


ISO 


I 


Integrating the above equation we get, where for whole region 

= O , 


2 ' 2 ^ 

i Cw+1)3 -% (2-0 y r (r-i) q/rK w , a,6> 


fc =5 r 


1 J. 


Where 


= He/ K and K is constant of integration. Wi th the 


aid of equation <4> , we can wright . 


U T+l ^i 2 {w-l> + 0 

= Ci+ Kr 

K ^ 


•7 2 w 

+ 0^ (w+1) -1(2-0 y ir-i) r q/Tk w 

e. 1 


<i7> 


w 


U 


— = Ci+ Kr 

a 4 

a 


T+i |,C2iw— i)+ (2 si) i T — 1) r qy k wt 


(18) 


^ t r" 1 @ 1 Cl LISS O "f SjIlOCl*-. 

Sililiiariy for a strong magnetic tiei ^ . 

CQii'iiitiofi C 5 ) iri ( 7 ) we get 


fP 


d€ 

■, e 


'( i- 


o ■ T , 

dP ' 2 

2 o 

2H H . 
o ' o 


dH ■ dfe 

p o d ^ 

^ b ,, 

■ ./G'... 7.. Q 


q dr 


H 


o 
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ir-iy 


dr 


<19> 


H 


dH / H and db / b in the equation 


Substituting the value of dP , „ „ 

o .. a Q Q O,, 

(19) and integrating we gat 


<3w— 1) 

- j~ 

k r 0^- 




+ 2(l-w> 3 0 ) (T“l>/wk 0 


< 20 ) 


where k is cons 
write 


tant of intqratian.mth the equation of (5) we 


(3w-i) 

U 3 - -X“ ^ 

— =Ci+ kr 0 

K 4 


w 


2 


3 — ibCi"' ^ 0 >qr (T^ll/wk 0 T 
*e ^ 


( 21 ) 


U 


= Ei-i- Kr 

4 


w 2 2 

.cifr:” +2 )qr i-r-l>/»k 

0^ *e 

( 22 ) 


-or strong shock ^if we 
E7? we get " 


substitute the shock condition (6) m 


dP 


r — ^ ( — ) 3 dU +E 

^ ■ H T P 


o 


5 , da 

+ % 1 

^ a crt-D + iiit'-J? ^ 

o' 





Of 



ios&am 



>'^0 tw4“i ) C2-P ^ 

e , 


1 


Cw+i^ 1 ■ ^2 ii ,i 




irn^l} 1 -% q 


%E1— ^ !8 3 Qf” ('t-D/TWK. {3 


< 3 w-i j /P 


-2<1-3w) 3 ) ) qr tT-1 > /Twk Ji 


,i ii^+i> 3 -% >j q»' 


('T"l)/'TWfc 0 


( 28 ) 


shock 


— Cis| 4 'BI 





Hw/ m+Bi 


,, = ). (K; r 


D 'qr 


ha ahpraasion (17) and (21) rapradant tha prnpaQatlon of »aak 

l.varg^n, cylindrical ahoc. wava through a rotating gaa with 

• j. +-ho nrja^ence of transverse weak 

(ftect of artificial viscosity in the p 

,nd strong magnetic field respectively - 

In relation <25> first term represents the solution obtained 

the i 4 ithai«, showing the effect of density distribution . Ihe 

■ 1 =i-r-™int of the corioiis force and 

second term arises obeviousiy on acco . . 

tha masnatic fiald . Tha valocity of aagnatogaadynaa.c shock 
..va is dat.raihad by first tern only when aagnatic fiald is 
strong , for larga valnas of r , tha first tram ts very small in 


strong 


;omoarison to the second term 


As a consequence, the shock velocity attains a mi ' 

value for a certain propagation distance Payond which the cpriplis 

tram stand dpminating and shock propagation distance 

.ci.iaL t-i-f mpa shock velocity may be 
i - 4 -»n«i sniniBium vaiUe qt i-ne - 


trem stand dominating an? 


corresponding to the minimum vo 
obtained by equating 
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